INVERSE ESTIMATES FOR ELLIPTIC BOUNDARY INTEGRAL 
OPERATORS AND THEIR APPLICATION TO THE ADAPTIVE 
COUPLING OF FEM AND BEM 



M. AURADA, M. FEISCHL, T. FUHRER, M. KARKULIK, J. M. MELENK, AND D. PRAETORIUS 

Abstract. We prove inverse- type estimates for the four classical boundary integral oper- 
ators associated with the Laplace operator. These estimates are used to show convergence 
of an /i-adaptive algorithm for the coupling of a finite element method with a boundary 
element method which is driven by a weighted residual error estimator. 



1. Introduction 

A posteriori error estimation and adaptivity have a long history in finite element meth- 
ods (FEMs), which harks back at least to the late 1970s. While the early mathematical 
analysis (see, e.g., the monographs [AO, BS, V]) focussed on a posteriori error estimation, 
significant progress has been made in the last decade in the analysis of the adaptive FEM 
(AFEM) regarding provable convergence and achievable convergence rates. For linear ellip- 
tic model problems and discretizations by fixed order polynomials on shape-regular meshes, 
convergence and even quasi-optimal convergence rates for AFEM have been proved; we refer 
to [CKNS] for symmetric problems as well as to [CN] for nonsymmetric problems and the 
references therein. 

The situation is less developed for the adaptive boundary element method (ABEM) and 
even worse for the adaptive coupling of FEM and BEM. For the ABEM based on first kind 
integral equations, [FKMP] and [G] proved very recently convergence and optimality. Specif- 
ically, [FKMP, G] studied lowest order discretizations of equations related to the simple-layer 
operator and hypersingular operator. We highlight that the symmetry of these two operators 
is an important ingredient in the optimality proofs for the ABEM in [FKMP, G]. As a first 
step towards a full analysis of the more complex case of the adaptive coupling of FEM and 
BEM, we show in the present paper convergence for Costabel's symmetric coupling. 

Broadly speaking, the procedure in [FKMP, G] and the present work relies on a framework 
delineated for AFEM in [CKNS]. The starting point for AFEM are reliable residual type 
error estimators. Ideally, the residual is measured in a dual norm; for example, in the 
classical Laplace Dirichlet problem with numerical approximation one has to evaluate 
||/+Aitft||i3--i(n)- Since such duals norms are difficult to realize computationally, the classical 
residual error estimators mimic local versions of them by weighted L 2 -norms of various 
components of the residual. The appropriate weight is given in terms of the local mesh 
size function h. Returning to the example of the Laplace Dirichlet problem, these are the 
elementwise volume residuals \\h(f + Auh)\\L 2 (T) an d the edge/face jumps of the normal 
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derivative ||/i 1//2 [<9 n Wh] Hx^f jst) - 111 effect, the residual is measured in stronger, but /i-weighted 
Sobolev norms. Inverse estimates are therefore a key ingredient to showing efficiency of 
these estimators and convergence in the context of adaptive methods. A second feature of 
error estimators in the FEM is the local character of the volume and edge/face residuals. 
This feature stems from the fact that a differential equation is considered. As a result, the 
only inverse estimates needed in the FEM are classical ones relating stronger integer order 
Sobolev norms to weaker ones. 

Our point of departure for ABEM are computable weighted residual error indicators made 
available in [C2, CMS, CMPS, CS1]. Analogously to the FEM, the nonlocal nature of the 
norm in which to measure the residual (these are typically non-integer Sobolev norms) is 
accounted for by /i-weighted integer order Sobolev norms; for example, ||/i 1//2 Vr(-)Hi, 2 (r) is 
taken as a proxy for | • |#i/2( r ), where Vr is the surface gradient on the surface T. The second 
feature of the residual error estimators in FEM mentioned above is the local character of 
the residual. This feature is not present in boundary integral equations. For example, the 
weighted residual error estimator for Symm's integral equation involving the simple-layer 
operator V is ||/i 1//2 Vr(/ — Vuh)\\L 2 (r)', even when restricting the integral to a single element, 
the nonlocal nature of V involves the numerical approximation Uh on whole surface T. As a 
result, the classical inverse estimates for spaces of piecewise polynomials, which were suitable 
in the FEM, are insufficient for the BEM. The appropriate inverse estimates are provided in 
the present paper. In this connection, the works [FKMP, G] are particularly relevant. Using 
similar techniques, [FKMP] considers the special case of the lowest order discretization of the 
simple-layer operator QJ, whereas the present article covers arbitrary (fixed) order conforming 
discretizations of all four operators. The work of [G] leads to very similar results; possibly 
due to the use of wavelet techniques, [G] requires C 1,1 -surfaces. Instead, our analysis relies on 
techniques from local elliptic regularity theory, and we may thus admit polyhedral surfaces 
here. We finally note that all four operators appear in our formulation of the FEM-BEM 
coupling. 

The remainder of this work is organized as follows: Section 2 collects all necessary notations 
and preliminaries (Section 2.1-2.2) and proves the new inverse estimates (Theorem 1 and 
Corollary 2), which are the main results of this work. Our analysis relies on elliptic regularity 
estimates for the simple-layer potential 03 (Section 2.4) and the double-layer potential ^ 
(Section 2.5). In Section 3, we consider an adaptive algorithm for Costabel's symmetric 
FEM-BEM coupling, which is applied to a linear transmission problem (Section 3.1). Our 
discretization includes the approximation of the given data so that an implementation has 
to deal with discrete boundary integral operators only (Section 3.4). We therefore extend 
the reliable error estimator of [CS2] to include data approximation terms in Proposition 13. 
Adapting the concept of estimator reduction [AFLP], which has also been used for (h — h/2)- 
type error estimators in [AFP], we prove that the usual adaptive coupling (Algorithm 14) 
leads to a perturbed contraction for the error estimator ge and thus obtain convergence ge — > 
as i — > oo. Since ge provides an upper bound for the Galerkin error, which unlike [AFP] 
does not rely on any saturation assumption, we thus obtain convergence of the adaptive 
FEM-BEM coupling (Theorem 15). A short Section 3.8 discusses the extension of these 
result to nonlinear transmission problems. Numerical experiments in Section 4 illustrate the 
convergence of the adaptive FEM-BEM coupling procedure and give empirical evidence that 
the optimal order of convergence is, in fact, achieved. 
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2. Inverse estimates for integral operators 

In this section, we prove certain inverse estimates for the four classical boundary integral 
operators associated with the Laplace operator. Independently of our work, similar estimates 
for fixed order piecewise polynomials have been proved in [G] by means of wavelet-based 
techniques. For technical reasons, [G] assumes the boundary to be fairly smooth, namely, 
C 1 ' 1 . In contrast, the present analysis is based on PDE techniques and allows us to treat 
polygonal/polyhedral boundaries. We mention that the estimate (2.12) below has already 
been shown in our earlier work [FKMP] in a discrete setting for lowest order elements £ 
V°(£e) and is generalized here to the case of arbitrary ip £ L 2 (T). 



2.1. Preliminaries & general assumptions. Let n be a bounded Lipschitz domain in 
M. d , for d > 2, with polygonal/polyhedral Lipschitz boundary T = dQ. The exterior unit 
normal vector field is denoted by v. We define f2 ext := IR d \ Q. Let Si denote a conforming 
triangulation of T into simplices, i.e. 

• Si is a finite set of non-degenerate (d — l)-dimensional compact surface simplices, i.e., 
affine images of the reference simplex E Te f C with positive surface measure; 

• r = LW £ E > Le -' F is covered by S e ; 

• for each E, E' £ Si with E ^ E', the intersection E l~l E' is either empty, or a j- 
dimensional simplex for j — 0, . . . , d — 2, i.e., a joint node, or a joint edge, or a joint 
face, etc. 

Moreover, we assume that Si is K-shape regular, i.e., 

diamfE)^ 1 n , . , 

max ^ < k < oo for d > 3 (2.1) 

Ee£ e \E\ ~ y ' 

with | • | the {d— l)-dimensional surface measure and diam(-) the Euclidean diameter, whereas 

max diam (^) < K < oo for d = 2. (2.2) 
e,e'gSi diam(E') 

Note that, for <i > 3, conformity of ^ and K-shape regularity (2.1) also imply (2.2) (though 
with a different, but bounded constant k). With each triangulation S#, we associate the local 
mesh size function hi £ L°°(r) which is defined elementwise by Hi\e '■= hi(E) := \E\ 1 '( d ~ 1 ^ 
for all E £ Si. We stress that K-shape regularity of Si implies Hi\e — dia,m(E). 
Let 7 C T be a relatively open subset of T. With 

Sj := S{\~{ := {E £ Si : EC. 7}, (2.3) 

we denote the restriction of Si to 7. It is always assumed that 7 is resolved by Si, i.e., Sj is 
a /t-shape regular and conforming triangulation of 7. 



2.2. Sobolev spaces and boundary integral operators. In this section, we very briefly 
fix our notation concerning Sobolev spaces and boundary integral operators and refer the 
reader to the monographs [M, HW, SS, Ver84] for further details and the precise definitions. 

For the boundary V = dQ of the bounded Lipschitz domain Q, we denote by Vr(-) the 
surface gradient. The Sobolev space i? 1 (r) can be defined as the completion of the Lipschitz 
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continuous functions on T with respect to the norm 

IMItfi(r) := IMlW) + l|Vr«||ia(r)- 

We denote by 7o nt (-) the interior trace operator, i.e., the Y^u is the restriction of a function 
u £ to the boundary T. The space H l / 2 (Y) is the trace space of H l (VL) equipped 

in the standard way with the quotient norm. For relatively open subsets 7 C V and s £ 
{-1/2,0,1/2}, we denote by 

#i/2+«( 7 ) = | w | 7 . v e frV2+-(r) jB uppi; c 7}. 

the space of all restrictions of functions to 7 and endow this spaces with the corresponding 
quotient norms. In particular, if v £ H 1 ^ 2+s {^) is extended by zero to the entire boundary 
r, then v £ H^+'iT) and || v ||^i/2+s^) — \\v Hij-i/a+sfr)- Finally, negative order spaces 

H' 1 ' 2 ^) := H^ 2 (T)', H~^ 2+s \-f) := # 1/2+s (7)', and H^ 1 ' 2 ^^) := H 1 ' 2 ^^)' 

are defined by duality, where duality pairings (-,-} are understood to extend the standard 
L 2 -scalar product. We note the continuous inclusions 

H ± d/2+-)( 7 )c J ff ± (V2+-)( 7 ) as well as h^ 1 / 2+s \Y) = H ±{ - l ' 2+s \Y). 
The interior conormal derivative operator 7] nt : H\{Vt) — > H~ l l 2 (T), where 

Hi(n) := {u £ H\Q) I - Am £ L 2 (fi)}, 
is defined by the first Green's formula, viz., 

( 7 l nt n , v) r = (Vu , Vu) - (-An , u) n for all v £ tf 1 ^). (2.4) 



Remark. The operator 7™* generalizes the classical normal derivative operator: if u £ 
H\{VL) is sufficiently smooth near a boundary point xq, then 7} nt can fre represented near x 
&?/ a function given by the pointwise defined normal derivative d u u. □ 

The exterior trace 7g xt and the exterior conormal derivative operator 7j xt are defined 
analogously to their interior counterparts. To that end, we fix a bounded Lipschitz domain 
U cR d with ncU. The exterior trace operator % xt : H X {U \ O) i/ 1 / 2 (r) is defined by 
restricting to T, and the exterior conormal derivative 7^ xt is characterized by (jf^u , v)r = 

For a function tt that admits both derivatives, we define the jumps [jiu] := 7^ xt u — Y^u, 
as well as [u] = 7§ xt u — 7Q nt u. 

We denote by G the fundamental solution of the d- dimensional Laplacian 

I — t^tt log la; — 2/1, for d = 2, 
G(x,y) = < |s ' 2.5 

\ + ^|x-y|-( rf - 2 ), forrf>3, 1 J 

where IS 01-1 ! denotes the surface measure of the Euclidean sphere in M. d , e.g., |§ x | = 2tt and 

An. The classical simple-layer potential 53 and the double layer potential & are defined 
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by 

Q3if>)(x) := J G(x,y)i>(y)dT(y), (kv)(x) := J d v{y) G(x,y)4>(y)dr(y), x G R d \ T; 

here, d u ( y ) denotes the (outer) normal derivative with respect to the variable y. These 
pointwise defined operators can be extended to bounded linear operator with 

^eL{H- ll2 {T);H\TJ)) and R G L(H 1/2 (T); H\U \ V)) . (2.6) 

It is classical that M5ijj = = AMv in U \ T for all ip G H- 1 ' 2 {Y) and v G H 1/2 (T). The 
simple-layer, double-layer, adjoint double-layer, and the hypersingular operators are defined 
as follows: 

QJ = 7 i nt 2T, £ = _ + # = -I + 7 i n %, and 2H = -7^. (2.7) 

As is shown in [HW, M, Ver84] , these linear operators have the following mapping properties 
for s G { — 1/2, 0, 1/2} and representations (in the case of the hypersingular operator 2U, the 
integral is understood as a part finite integral): 



SJG L(^ 1 / 2+s ( 7 );i/ 1 / 2+s (7)), 


(QJ^)(x) 


= j G(x 7 y)^(y) dT(y), 


(2.8) 


^eL(ff 1/2+8 (7);^ 2+! ( 7 )), 


(^«)(z) 


= J d„(y)G(x,y)v(y) dT(y), 


(2.9) 


# G ^(if- 1 / 2 ^^);^- 1 / 2 ^^)), 


(£ty>)(a:) 


= J d u{x) G(x, y) v(y) dT(y), 


(2.10) 


OH G L(i/ 1 / 2+s ( 7 ); _fT~ 1 / 2+s ( 7 )), 


(22Jt;)(x) 


= -du(x) / d u{y )G(x,y)v(y)dT(y). 


(2.11) 



Jr 



2.3. Statement of main result on inverse estimates. The following theorem is the 
main result of this work and the mathematical core of the arguments that allow to transfer 
convergence results from AFEM to ABEM. 

Theorem 1. There exists a constant C\ m > such that the following estimates hold: 



/4 /2 V r 2¥lb( 7) < CiBvOHItf-v^) + ll^Vlbw), (2-12) 

\\h) /2 ^\\L^) < Cinv(||^b- 1/2(7) + H^Vll^w), (2.13) 

\\h\ /2 V v &v\\ L , h) < C inv (\\v\\ ni/2{1) + \\h l /\ r v\\ L2M ), (2.14) 

\\h] /2 mv\\ L 2 {l) < CW(IMIhV*( 7 ) + ||^ 1/2 Vr^||L 2 ( 7) ), (2.15) 



for all functions ip G L 2 ^) and all v G H 1 ^). The constant C inv > depends only on T, 7, 
and K-shape regularity of £g. 

Remark. The estimates (2.12)-(2.15) are rather easy to show for globally quasi-uniform 
meshes, i.e., hg{E) ~ hi(E') for all E, E' G £i. For example, to see (2.12), one recalls 
stability of ^3 : L 2 ( 7 ) H 1 ^) to get 

||^ 1/2 Vr^l|L 2(7 ) < ^ 1/2 ||^||^(7) £ h\ /2 U\\m,) * 11^11^(7)- (2-16) 
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One sees that (2.16) is slightly stronger than (2.12), where an additional term llV'lljf-i/a^) 
arises on the right-hand side. □ 

For each element E £ £J, let '■ E m f — > E denote the affine bijection from the reference 
simplex E Tef C IR^" 1 onto E. We introduce the space of (discontinuous) piecewise polynomials 
of degree q and the spaces of continuous piecewise polynomials of degree p by 

V q {£j) := {^e £ L 2 {~f) : V£ £ £j ^ o lE is a polynomial of degree < q}, (2.17) 

: = V p (£]) n C( T ), and S P (£7) := {V<| 7 : ^ £ 5 P (^) with supp(^) C 7}. 

(2.18) 

We note the inclusions V q [£j) C L 2 ( 7 ) C ^~ 1/2 (7), Sl{£j) C tf 1 ^) C # 1/2 (7), and 
S p (£j) C H 1 ^), as well as So(£e) = <S p (£e) in case of 7 = T. If we now restrict the 
estimates (2.12)-(2.15) of Theorem 1 to discrete functions £ V q {£j) and £ Sq(£J), we 
obtain the following estimates. 

Corollary 2. There exists a constant C- mv > such that the following estimates hold: 



max{l,g}-l/4 /2 V r 2I^|U2 (7) < Ci„ v ||^|| 5 _ 1/2(7) , (2.19) 

max{l,g}- 1 ||/ i J /2 ^|| x2(7) < C inv ||^|| s _ 1/2(7) , (2.20) 

p-^hyWrmh^y) < C- mv \\Ve\\ Sl/2{jV (2.21) 

p- 1 ||/ i J /2 2UT4|| i 2 (7) < a nv ||^|| #1/2(7) , (2.22) 



/or a// discrete functions £ V q {£j) and Vi £ Sq(£]). The constant C- irw > depends only 
on 7 and i/ie shape regularity constant k of £e, but is independent of the polynomial degrees 
q > resp. p > 1. 

Remark. We stress that for discrete spaces with locally varying polynomial degrees, e.g. 

V qi {£j) := {^e £ £ 2 (t) : V£ £ £/ ^ o^ E is a polynomial of degree < qi(E)}, 

the inverse estimates of Corollary 2 remain true as long as the polynomial degrees are compa- 
rable on neighboring elements, i.e., qe{E) ~ qe{E') for E R E' 7^ 0. For example, the inverse 
estimate (2.19) involving QJ then reads 

||/i; /2 max{l,g £ } _1 Vr^|U 2 ( 7 ) < a nv ||^|| 5 _ 1/2(T) for all £ V q *(£]), 

and the estimates (2.20) -(2.22) can be extended analogously. We refer the reader to [K, 
Theorem 4.4], where the inverse estimates are proved in this extended fashion. □ 

Proof of Corollary 2. Our starting point are two inverse estimates from [EG, Theorem 3.9] 
and [AKP, Proposition 3]: 

max{l, g }- 1 ||^ 1/2 ^|| L2(7) < ||^|| g _ 1/2(7) for all ^ £ V q {£]), (2.23) 

p- 1 ||/i] /2 V r ^|| i2(7) < |N| 5l/2(7) for all V, £ S p (£]), (2.24) 

where the hidden constants depend solely on T and the K-shape regularity of £J. Combining 
(2.23) with (2.12)-(2.13) leads to (2.19)-(2.20); the bound (2.24) in conjunction with (2.14)- 
(2.15) yields (2.21)-(2.22). □ 
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The remainder of this section is devoted to the proof of Theorem 1. The proof will first 
be given for 7 = T, and the general case 7 C T is deduced from it afterwards. 

On the technical side, an important difficulty of the proof of Theorem 1 arises from the 
fact that the boundary integral operators are nonlocal. We cope with this issue by splitting 
the operators into near field and far field contributions, each requiring different tools. The 
analysis of the near field part relies on local arguments and stability properties of the BIO. 
For the far field part, the key observation is that the BIOs are derived from two potentials, 
namely, the simple-layer potential QJ and the double-layer potential R by taking appropriate 
traces. Since these potentials solve elliptic equations, inverse type estimates ( "Caccioppoli 
inequalities") are available for them. The study of these two potentials is the topic of 
Sections 2.4 and 2.5. We will discuss the case of the simple-layer potential in greater 
detail first and be briefer afterwards in our treatment of the double-layer potential ^, since 
the basic arguments are similar to those for QJ. 

2.4. Far field and near field estimates for the simple layer potential. We start with 
subsection 2.4.1, where we introduce the decomposition of the simple layer potential into far 
field and near field. For either of this parts, we provide inverse estimates. Subsection 2.4.2 is 
devoted to the derivation of inverse estimates for the near field parts, whereas we deal with 
far field parts in subsection 2.4.3 

2.4.1. Notation and decomposition into near field and far field. For each element E G Se and 
5 > 0, we define the neighborhood Ue of E by 

Ed U E := |J B 2Shc(E) (x), (2.25) 

where B £ {x) := [y G W 1 \ \x — y\ < C IR d . By K-shape regularity of £1, there exist 5 > 
and MeN such that T D Ue is contained in the patch ue{E) of E, i.e., 

r n U E C u t (E) :=\J{E' eS e : E'nE^®}, (2.26) 

and that the covering T C {J E ee Ue is locally finite, i.e., 

#{U E : E G 8 t and x G U E } < M for all x G R d and £ G N. (2.27) 

Finally, we fix a bounded domain U <Z^ d such that 

U E C U for all E G S t . (2.28) 

To deal with the nonlocality of the integral operators, we define for functions ip G £ 2 (r) and 
E G Si the near field and the far field of the simple-layer potential = by 

:= %U> X mu B ) and u% E := ^ X r\u E ), (2.29) 

where Xu denotes the characteristic function of the set 00 C IR d . We have the obvious identity 

u<a = u ^ + u^ E for all E G S t . (2.30) 

In our analysis, we will treat wJjff and u^ r E separately, starting with the simpler case of 



2.4.2. Inverse estimates for the near field part uf^™. The near field parts of a potential can 
be treated with local arguments and the stability properties of the associated boundary 
integral operators. 

Lemma 3. There exists a constant C near > depending only on T and the n-shape regularity 
of Ei such that for arbitrary E £ Ei and ^>f £ V°(£e) with supp (^f) C u>e(E) there holds 

E \\^f\\h(p B ) < CU r £ \\hl /3 9f\\h( Ut wy 

EeS e Ee£ e 

Proof. We fix an element E £ £(_. Denoting by ^ff(E') the value of \I/f £ V°(£e) on the 
element E' we compute 

(Vafaf = E *F(^') / V s Gf(s, y) dT(y) for all x £ R d \ T. 

The number of elements E' in the patch u^{E) is bounded in terms of the shape regularity 
constant k and thus 

|(w*f)(*)| 2 < E l*f(^)l 2 (/ ; |v,G(x,y)|dr(y)) 2 

with some E'-independent constant, which depends only on k. Since the mesh £% is re- 
shape regular, we can select a constant c > 0, which depends solely on k, such that UgC 
B C h e (E)(bE') and £" C B che ( E ){bE>) for each neighbor £" £ ui(E). We integrate over [/# and 
estimate the remaining integral 



/ |(VQJtff)0r)| 2 ^< E l^')| 2 / (/ |V,G(x,y)|rfr(y)) 2 dx 

< E i*iw/ (/ —^r( y ))V 

E'^ t {E) J B chdE )(b E >) JB chf(E) (b E ,)nr E , P </l 



where T^/ denotes the hyperplane that is spanned by E' . Scaling arguments then yield 
/ \{VWf){x)\* dx < E \*f(E')\ 2 h e (E) d [ (f | ^dr^Vd* 

JI7k c,^ . J Bi(0) v >/_Bi(0)nR d - 1 P — 2/| 7 



< e i*fwrw l ^ii^ 1/a « r Fii^( W4 ( B ))- 

Summing this last estimate over all E £ £g, we conclude the proof. □ 

Proposition 4 (Near field bound for 93). There is a constant C ncar > depending only on 
T and the n-shape regularity of ' £(_ such that the near field part u^ 3 ^ satisfies u^ 3 ^ £ H l (U) 
and 7i nt M^E £ H X {T) as well as 

E ll^v^S 2 -^ + E \^ u Te\\Ihu e) < cu ll^Vlli-cn- ( 2 - 31 ) 

Ee£e Ee£ e 
Proof. The stability assertion 93 : L 2 (T) — y H l (T) proved in [Ver84] gives, for each E £ £e, 
HVr^K^IU^) < \\n^Xu E nr)\\m { r) < Uxu E nv\\^(r) = IMU^nr)- 
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Summing the last estimate over all E G £g and using the finite overlap property (2.27) of 
the set Ue, we arrive at 

ll^ /2 Vr7i nt ^ll^(^) = E ^(^)llVr7i nt ^lli 2(£ ) 

Eee e Ees e 

<E^^)^HW B nr)^llOlli 2( r), 

EG£ e 

where all estimates depend only on the K-shape regularity constant n. This bounds the first 
term on the left-hand side of (2.31). To bound the second term, let denote the L 2 (T)- 
orthogonal projection onto V°(£e). We decompose the near field as u% a E = QJ(n^(^Xrn(7 B )) + 
QJ((1 - Il e )ipxmu E ) ■ The condition supp(^Xmc/J C u e (E) implies supp (U e (^xmu E )) C 
ue(E) and therefore, taking tyf = Uf(if)Xrnu E ) i n Lemma 3, we conclude 

E !|V^(n,^Xrnt/J)||i 2(c/B) < E \\h l / 2 M^Xrn Uls )\\% M E)) % ll^fVll W), (2-32) 

Ee£ e Ee£ e 

where we used the local L 2 -stability of Hi in the last estimate. Recalling the stability 
«U : H~ l ' 2 {Y) H l {U) of (2.6) we can estimate 

|| W((l - n £ )^ X rnt/ B ) \\^u B ) < || W((l - Ii^xmu E ) \\l^u) 

< ||(1 - lWxrnaJlif-i/2(r) 
Together with a local approximation result for IL? from [CP1, Theorem 4.1], we obtain 

-n e )i} X rnu E )\\h {UE) < E W^-^XmuJl-^ 

Ee£ e E££ e 

^ E ll^ /2 (V ; Xrnc/ E )|li2 (r) ~ ||^ 1/2 ^|li2 (r) . (2.33) 
Ee£ e 

The combination of (2.32)-(2.33) yields the desired estimate in (2.31) for ^2 E€£( \\^ u ^e\\l 2 (u e )- 

□ 

2.4.3. Estimates for the far field part u^ r E . The following lemma is taken from [FKMP]. 
For the convenience of the reader and since the same argument underlies the proof of the 
analogous lemma regarding the double-layer potential, we recall its proof here. 

Lemma 5 (Caccioppoli inequality for u%j T E ). There is a constant C cacc > depending only 
on the n-shape regularity of Si such that for the function u^ E of (2.29) the following is true: 
u&ln e C°°{Q), ugfrsln- G C°°(fi cxt ), and u^ e \ Ue G C°°(U e ) with 

\\D 2 u%] E \\ L 2 {BshdE){x)) < C cacc — — \\Vu^ E \\ L 2 [B2Shmix)) for allxeEe S e . (2.34) 

Proof. The statements u^ E \ n G C°°(fi) and w|^|a-* G C°°(fi cxt ) are taken from [SS, Theo- 
rem 3.1.1], and we therefore focus on the statements on u^ e \ Ee G C°°(U e ) and the estimate 
(2.34). According to [SS, Proposition 3.1.7], [SS, Theorem 3.1.16], and [SS, Theorem 3.3.1], 
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the function u% T E G H} oc (R d ) := {v : R d -> R : v\ K G /^(if) for all A" C M d compact} 
solves the transmission problem 

a.e. in Q U fi cxt 

in iyV2(r) (2.35) 

-^Xr\u E inff-V2(r). 

In particular, (2.35) states that the jump of wffe as well as the jump of the normal de- 
rivative vanish on r D U B - This implies that uyj,E is harmonic in Ue by the following 
classical argument: First, we observe that u%j T E is distributionally harmonic in Ue, since a 
two-fold integration by parts that uses these jump conditions shows for v G C^°(Ue) that 
(u%f E ,—Av) = 0. Weyl's lemma (see, e.g., [MO, Theorem 2.3.1]) then implies that u% r E is 
therefore strongly harmonic and u% E G C°°(Ue). 

The Caccioppoli inequality (2.34) now expresses interior regularity for elliptic problems. 
Indeed, [MO, Lemma 5.7.1] shows 

\\D 2 u\\ L 2 {Br) < (\\f\\ L 2 (Br+h) + - \\Vu\\ L 2 (Br+h) + ^\\ u \\ L 2 {Br+h ^ (2.36) 

for each u G H l (B r+ h) such that u G H 2 (B r ) and Am = / on B r+ h with balls B r C B r+ h 
with radii < r < r + h and some / G L 2 (B r+ h); the hidden constant depends solely on the 
spatial dimension and is independent of r, h > 0, and u, f. We apply (2.36) with / = and 
M = u k,e ~ c e, where c s = ^ f B26h (B) ( X ) Usin g additionally a Poincare 

inequality then leads to (2.34). □ 

The nonlocal character of the operator QJ is represented by the far field part. Lemma 5 
allows us to show a local inverse estimate for the far field part of the simple-layer operator: 

Lemma 6 (Local far field bound for QJ). For all E G £g, there holds 

||/4 /2 Vr7^<*IU*( B) < ll^ /2 V< s |U*( fi ) < C far \\Vu^e\\lhu e ). (2.37) 
The constant Cf ar > depends only on F and the K-shape regularity constant of Si. 

Proof. By Lemma 5 we have u l ^ E G C°°(Ue). The first estimate in (2.37) follows from the 
fact that, for smooth functions, the surface gradient Vr(-) is the orthogonal projection of 
the gradient V(-) onto the tangent plane, i.e., Vr7o nt w(a0 = Vu(x) — (V«(x) • v(xj) u(x), 
see [Ver84]. 

To prove the second estimate in (2.37), we fix an E G Si. First, we select N points Xj G E, 
j = 1, . . . , N, such that 

N 

EC [jBsh^ixj). 

j'=i 

We may assume that the number N of points depends solely on the K-shape regularity of 
the mesh. This follows by geometric considerations as detailed in [FKMP, Lemma 3.5]; 
essentially, one may select the points Xj from a regular grid with spacing ^5hi(E) and cover 
E with balls of radii 5hz(E) centered at these points. The centers outside E that are required 
for the covering are then projected into E to ensure that all points are in E. 
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hi 



far 1 



Next, let Bi := Bsh t {E){xi) and Bi := B 2 sh e (E){xi) C U E - Using a standard trace inequality 
and the Caccioppoli inequality (2.34) we infer for all indices % that 

f 2 1 f 2 f 2 f 1 f 2 

HVUjijyiia^nE) < ^T^yllVw^ll^^ + || Vw^ B || L 2 (B .) ||D m^|| L 2 (B]) < ——. [^1^11^^. 
We use the last estimate to get 

N 1 JV 

II V7 , int „ . f a.r 1 1 2 ^ ||Y7„,far II 2 s~ \ ^ II Y7„,f ar II 2 < \ ll\~7„,f ar II 2 

H v r7o U E IIl 2 (E) ^ II V u <S,e\\L 2 (E) ^ ^ II V %J,£llL2(B;n£) ^ T / , || vm QJ,eII/-,2(-r,) 

i=l ^ ' i=l 

This concludes the proof of (2.37). □ 
Summation of the elementwise estimates of Lemma 6 yields the following result: 

Proposition 7 (Far field bound for 23). There is a constant Cf ar > depending only on T 
and the K-shape regularity of Se such that 

\\hl /2 V rl ^u%Jl 2{E) < ll^ /2 V< E ||i 2(£) < Q ar (Wi'Wl-^r) + WK^WUd] 
Ees e Ees e 



Proof. We use the local far field bound (2.37) of Lemma 6 and u^ E = %Jip — u% E ', 

\l 2 (Ue) 



The first term on the right-hand side in (2.38) is estimated by stability of 03 and the finite 
overlap property (2.27) 

The second term in (2.38) is bounded with the aid of the near field bound (2.31). □ 



2.5. Far field and near field estimates for the double layer potential. Section 2.4 
studied the simple layer potential in detail. Corresponding results for the double layer 
potential are derived in the present section. 

2.5.1. Decomposition into near field and far field. We use the notation introduced in Sec- 
tion 2.4.1. Additionally, in order to define the near field and far field parts for the double-layer 
potential, we need an appropriate cut-off function: Let fife denote the set of nodes of Se- For 
any E G Se define 

z&M e noj e (E) 

where r\ z G S 1 (Se) denotes the hat function associated with the boundary node z, i.e., 
r] z G S l (Se) is characterized by the condition r) z (z') = 5 zz i for all z' G fife, where 5 ZZ > 
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denotes the Kronecker delta. We use the abbreviation Qt(E) := ui(uJi(E)) := \J {E' G S £ : 
E' fl uji(E) 7^ 0} for the second order patch, where we recall from (2.26) that coe(E) denotes 
the patch of E G £(,. Note that u)i(E) = supp(r] E ) for the cut-off function r\ E of (2.39). We 
note 

VeIthUe = 1) VE\r\w t (u t (E)) = 0, ||^|U<x>(r) = 1, and ||Vrf7#IU°°(r) - h^E)' 1 , 

(2.40) 

where the constant involved in the last estimate depends only on the K-shape regularity of 
Si. For the double-layer potential u& = &v of a density v G H 1 ^) we define the near field 
and the far field part by 

u n £%:=&{(v-c E ) VE ) and u% E := I((v - c E ){l - r, E )), (2.41) 

where c E G K is a constant that will be specified below. Since 81 = — 1 in Q and 81 = in 
fl ext , we have, for every E G Si, the identities 

Ua + c E = v%% + v%* E infi and u st = v%$ + u% E in fi cxt . (2.42) 

2.5.2. Inverse estimates for the near field part u^ E . The proof of the near field bound for 
the double-layer potential needs an appropriate choice of the constants c E elin (2.41). 

Lemma 8 (Poincare inequality on patches). For given w G H l {T) and E G Si, there is a 
constant c E G R such that 

\\w - c E \\ L 2 {Qj , {E)) < C , i||^V r w||L2(S),(£;)), (2.43) 

1 /2 

\\(w - c jE )?7 jB ||^i/2 (r) < Ci||/i/ V r w\\ L 2 {Qe{E)) , (2.44) 

\\{w - c E )rj E \\ H i {r) < Ci||Vrw||i,2(D,(B))- (2.45) 

TTie constant C\ > depends only on the K-shape regularity constant of St and on the surface 
measure of T. 

Proof. The first estimate (2.43) is established by assembling local Poincare inequalities on 
u)i(E) with the help of [DS, Theorem 7.1]. The properties of the cut-off function r\ E detailed 
in (2.40), the estimate (2.43), and the product rule yield 

||V r ((w - c e )tie)\\u»(p) < - c E )VrVE\\L^r) + \\VEV r {w - c E )\\ L 2p) < ||VrHU 2 

since uit(E) = supp(r] E ). Hence, we obtain with the trivial bound hi(E) < jr) 1 /^- 1 ) < 1 

\\{w - c E )r) E \\ 2 H i {r) = \\{w - c E )r) E \\ 2 L 2 {r) + ||V r ((w - c E )r] E ) ||£ 2(r) < \\^rw\\l HQe{E)) . 

This proves (2.45). It remains to verify (2.44). To that end, we recall the interpolation 
inequality || ■ ||^ 1/2(r) < || • \\l^t)\\ • \\h^(v) and note that \\(w -c E )rj E \\ L 2 {T) < \\w - c E \\ L 2 {Se{E)) 
to get 

- c E )rj E \\ H i/2 (r) <\\{w- c E )rj E || ^ 2 2 (r) || (w - ce)^^^ 
< \\hiV r w \\% 2 mE)) || V r w \\^ Qtm 
- ll^ /2 VrHU 2 (£ f (£))> 

where the last estimate hinges on K-shape regularity of Si. □ 
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The following lemma provides an estimate for the near field part of the double-layer 
potential. 

Proposition 9 (Near field bound for £). Let v G H 1 ^) and consider u^ E defined by (2.41) 
with the constant ce given by Lemma 8. Then 7o nt w| e |; r G H 1 ^), u 1 ^'^^ G H l {Vt), and 
u^\u\Ti e H\U \Ti) with 

E (\\h\' 2 ^v^u^\\% {E) + llVti^ll^nn, + llV^Hi,^^)) 

E€£ t (2.46) 

The constant C ncar > depends only on V and the n-shape regularity of Si. 

Proof. First, the trace of the double-layer potential 7g nt ^ = & — | : i/ 1 (T) — > H 1 ^) is 
continuous, (2.9). Taking into account (2.40) and the Poincare-type estimate (2.45), we 
observe 

HVr7r« r |U 2( £) < ||Vr7r« r |U 2( r) < \\(v - c E ) VE \\^ tim) < \\ V r v\\ LHQe{E)) . 
Summation over all E G shows 

E ll^ 1/2 Vr7o nt « r |li 2(S ) < ll^ /2 Vr^||l 2(r) . (2.47) 

EeS/t 

Second, we use continuity of 8. : H l / 2 (T) -> H\U \ V) of (2.6) and get 

l|V^ r ||i 2([/B n^) + \\^<Ml^u E ^) % - c e )ve\\ 2 h i /Ht) < \\hl /2 V r v\\l 2{Qm) , 
where we have used (2.44) in the last step. Summation over all E G £(, gives 

E (llVti^lli-^nn) + llVu^Ui^no-)) < ll^ /2 V r .;||i 2(r) . (2.48) 
EeSt 

Combining (2.47)-(2.48), we conclude the proof. □ 

2.5.3. Estimates for the far field part u% E . As for the simple-layer potential, we have a 
Caccioppoli inequality for the double-layer potential, which underlies the analysis of the far 
field contribution. 

Lemma 10 (Caccioppoli inequality for u { ^ r E ). For the constant C cacc of Lemma 5 the func- 
tions u^ E of (2.41) satisfy v%* E \ n G C°°(fi) ; u| r E |n-t G C°°(fi ext ), and uf >E \u E G C°°(U e ) 
with 

\\D 2 uf iE \\ L 2 {BsH{E){x)) < C cacc — — \\Vu% E \\ L 2 {B2Shi{E){x)) for allxeEe S e . (2.49) 

Proof. The proof is very similar to that of Lemma 5. One observes that the far field u^ E 
solves the transmission problem 

-Auf E = a.e. in Q U fi cxt 

Wae\ = (v-c E )(l-r) E ) intfV^r) 
hxujfe] = infT-V2(r). 

We note that (1 — r) E )\mu B = by construction of r\ E in (2.39). Hence, the same reasoning 
as in the proof of Lemma 5 can be done to reach the conclusion (2.49). □ 
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Lemma 11 (Local far field bound for £). For all E e Eg there holds 

\\h) l2 V Tl ^u^ E \\ L 2 {E) < \\h\ /2 Vu^ E \\ LHE) < C far \\Vu^ e \\ l , {Ue) . (2.50) 

The constant Cf ar > depends only on V and the K-shape regularity constant of Eg. 

Proof. The lemma is shown in exactly the same way as the corresponding bound for the 
simple layer potential 03 in Lemma 6 appealing to the Caccioppoli inequality (2.49) instead 
of (2.34). □ 

Proposition 12 (Far field bound for There holds 

E II^VrT^II^) < E W^^eWUe) < Cfa r (||/4 /2 VHH 2(r) + IMI^ /2(r) 
Ee£t Ee£e 

The constant Cf ar > depends only on T and the K-shape regularity constant of En. 
Proof. We have by Lemma 11 

E \\ h l /2 vrir<E\\i HE) < E ii^v^iii^) < e \\^<e\\i H u b) 

E£S e Ee£z Ee£ e 



E 1 1 1 1 ^([/^nn) + E ll^ u ^ll£ 2 (E/finfi ext )- 

Ee£ e Ee£ e 



(2.51) 



With the properties in (2.42) and a triangle inequality we get 



E ll^ /2 V< E ||£ 2(B) < E (l|VS(« - c fi )|| 2 L2([7Bnn) + ||VJfe| 



|2 

lL 2 ([/ B nn ext ) 



4- f HT79/ near ll 2 -I- ll^77/ near ll 2 

+ yl VM £,£ \\L 2 (u E nn) + II v %,s llL 2 ([/ B nn cxt ) 



(2.52) 



The near field contribution is bounded by Proposition 9. Furthermore, noting VM.ce = 
V(— eg) = in f2, we get 

E ll^ 1/2v <sll^(^) S E (ll v ^lli 2 (c/ E nn) + liv^lll.^nnext)) + ||/4 /2 vHii 2(r) 

II 112 , II I,l/2v7 l|2 

z IpIIhi/2(d + II V v H| i2(r) , 

where we have used continuity of ^. □ 

2.6. Proof of Theorem 1 for 7 = T. We are now in position to prove the inverse 
estimates (2.12)-(2.15) of Theorem 1. 

Proof of the inverse estimate (2.12) for the simple-layer potential 2J and 7 = T. Let ip G L 2 (T). 
Then, 

H/^Vr^llW) = E II^Vr^Hi^) 

E££e , , (2-53) 

< E H^ 1/2 VrTo nt ^ll| 2 (i,) + E II^VrTo^^lllw 

Ee£ e EeSi 
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Both sums on the right-hand side can be estimated with the bounds of Propositions 4 and 7. 
This yields 

||/4 /2 V r 2^IU*(r) S IMI*-v at r) + llOlU^r). 
and concludes the proof. □ 

Proof of the inverse estimate (2.13) for the adjoint double-layer potential 8! and 7 = V. Let 
ip G L 2 {T). We split the left-hand side into near field and far field contributions to obtain 

\\hTm\W) < E HE)\W^Xu E nr)\\l> {E) + E MWtoWllW (2.54) 

The continuity # : L 2 (r) L 2 (r) stated in (2.10) yields for the near field contribution 
E h t( E )W&'^Xu E nr)\\h {E) < E HE)\\#(i>Xu E nT)\\h {r) < E M^ll^ll 

E££ e E££ e Ee£ e 

< ll^fVll^r)- 

For the far field contribution, we write u%j T E = %3(ipxr\u E ) an d no ^ e that 8! = — 1/2 + 7™*5J 
and clearly (iPx~t\u e )\e = 0. Therefore, on E we have A'(ipxr\u E ) = 7i nt ^ijV Furthermore, 
by the smoothness of u^ E near E (see Lemma 5) we have 

Ii'uke = duU^ E on E (cf. 

Remark 2.2) and get 

\\#Wxi\u B )\\i?{E) = WlT^WmE) = \\dvU% E \\ L 2 {E) < \\Vu^ E \\ LHE) . 
The far field contribution in (2.54) can therefore be bounded by Proposition 7 
E h e (E)\\&>^ XrxUE )\\i 2{E) < e !I^ /2 V< E ||| 2(S) < Wh^i'Wl^ + ii^ii^x 

Ee£i 

Altogether, this gives 



/2 (r) - 

EeSe EeSz 



||/iJ /2 i?V|| L2(r) < \\h 1 i /2 4j\\ L 2 {r) + ||^||h-i/2(d 
and concludes the proof. □ 

Proof of inverse estimate (2.14) for the double-layer potential 8. and 7 = V. Let v G iJ 1 (r). 
We recall the stability of R — \ + 7g nt jl : H 1 ^) — > H l (T), from which we conclude y^Siv G 
i/ 1 (r). Therefore, 

1 

2' 

(2.55) 



|/4 /2 Vr-MU 2( r) = ||^ 1/2 V r (| + 7i nt -ft)^||^(r) < - 1| V r ^|k-(r) + \\h\ ,2 Vvlt^\W) 



with = &v. There holds u& + eg = w^ e |; r + in f2, cf. (2.42). For the second term on 
the right-hand side in (2.55), we obtain with the constants ce of Lemma 8 that 

||/iJ /2 V r 7i n V||| 2(r) = E h e (E)\\V rl £\u A + c E )\\l HE) 

S ^ (2-56) 

< E h t {mVrT?«%B\\U*) + E Mmvr^u% E \\i 2{E) . 

E££ e EeSe 
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The first sum can be bounded by Proposition 9, whereas the second sum can be bounded 
by Proposition 12. Altogether, this yields 

||/iJ /2 V r ^||L2(r) < INI #1/3(1*) + II^VHIz^r) 
and concludes the proof. □ 

Proof of inverse estimate (2.15) for the hypersingular integral operator Q2J and 7 = T. Let a£ 
v G i/ 1 (r). We split the left-hand side of (2.15) into the sum over all elements E G En. On 
every element, we subtract the constants ce from Lemma 8. Note that Wc E = 0. Splitting 
now into near field and far field yields 

H^^llW) = J2 Whfmv ~ c E )\\h (E) 

< Whfmv - c e )ve)\\1 He) + £ W h i /2 W « v - c *)(l - VE))\\l HE y 
Ee£ e Ee£ e 

(2.57) 

The near field contribution is bounded by the stability of 2U : H 1 ^) — > L 2 (T) stated in 
(2.11): 

\\W((v - c E )r] E )\\ 2 L 2 {E) <\\(v- c E )r] E \\ 2 H1{Qem < \\Vrv\\ 2 L 2 {Qe{E)) , 

where we have used the Poincare-type estimate of Lemma 8 in the last step. The sum over 
all elements gives 

\\h l / 2 m(v - c E ) VE )\\l 2{E) < £ h e (E)\\V r v\\ 2 L2{Qem < \\h\ l2 V v v\\ 2 L2{r) . 

Ee£ e Ee£e 

It remains to bound the second term on the right-hand side in (2.57). In view of the support 
properties of r) E , the potential u^ E = &((v — c E )(l — r] E )) is smooth near E (cf. Lemma 10) 
so that 7i nt w| lI £ = d u Uft V E on E. Furthermore, since 20 = —7™*$. we see 

\m(v-c E )(l- VE ))\\ 2 LHE) = hi^&lll^) = \\d„u% E \\ 2 L2{E) < \\Vu% E \\ 2 L2{E y 
We use Proposition 12 to conclude 

\\hl /2 mv - c E ){\ - Ve))\\1 H e) < E ll^V^lli,^) < ll^ 1/2 V r ^||i 2(r) + INI^r)- 

Ee£ e E££ e 

Altogether, we obtain 

||/iJ /2 2U?;|| L 2 (r) < \\h\ /2 Vrv\\ L z(T) + IMI/fi/2(r) 
and thus conclude the proof. □ 

2.7. Proof of Theorem 1 for 7 ^ T. Finally, it remains to prove the inverse esti- 
mates (2.12)-(2.15) of Theorem 1 for the case 7 ^ V. Let ift G L 2 ^) and v G if 1 (7). We 
define the trivial extensions ip G L 2 (T) and v G H 1 ^) by 

~, x _ J r tp(x) if x G 7 — / \ f v(x) if x G 7 

W : ~ \ ifxGT\ 7 ' W W : " \ ifxGT\ 7 ' 
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Note that IMI 5-1/2 (7) = !MU-V2(r) and IMIi?V2( 7 ) = IMIhV2(t)- Witn tnis ; we see 

||/i; /2 Vr^||L2( 7) < ||/4 /2 V r 2¥||L2(r) < C™(U\\ H -i/2 ( r) + ||^ /2 * 2 (r)) 

= Cinv(||^ll#-i/2 (7) + Wh 1 / 2 ^]^), 

which proves estimate (2.12). The other estimates (2.13)-(2.15) follow with the same argu- 
ments. □ 

3. Convergent Adaptive Coupling of FEM and BEM 

In this section, we use the inverse estimates of Section 2 to prove convergence of an adaptive 
FEM-BEM coupling algorithm. 

3.1. Model problem. We consider the following linear interface problem 



-div(2lVu int ) 


= / 


in n int 


:= tt, 


(3.1a) 


-Au cxt 


= 


in fi cxt 


:= R d \Q, 


(3.1b) 


w int _ ^ext 


= Uq 


on r, 




(3.1c) 


(2lW nt - Vu cxt ) • v 


= 00 


on r, 




(3.1d) 


u cx \x) 


= 0(1/|S|) 


as \x\ - 


■> 00, 


(3.1e) 



see Section 3.2 below for some remarks on the radiation condition (3.1e) for 2D. Here, Q 
is a bounded Lipschitz domain in IR d , d = 2,3, with polygonal resp. polyhedral boundary 
T := dQ and exterior unit normal vector v. We assume that the symmetric coefficient matrix 
2l(x) G Mfy^ depends Lipschitz continuously on x and has positive and bounded smallest 
and largest eigenvalues 

< C min < X m i n (x) < A max (x) < C max < 00 for almost all x G Q (3.2) 

and x-independent constants C max > C m \ n > 0. The given data satisfy / G L 2 (Q), u G 
H l l 2 (T), and 0o G H~ l l 2 (T). As usual, (3.1) is understood in the weak sense, and the 
sought solutions satisfy u int G H\Q) and u cxt G H} oc (tt cxt ) = {v : fi cxt -»• R : v G 
H\K) for all K C Q cxt compact} with Vu cxt G L 2 {VL cxt ) d . 

With the boundary integral operators from (2.8)— (2.11), Problem (3.1) is equivalently 
stated via the symmetric FEM-BEM coupling, cf. e.g. [CS2, Theorem 1]: Find (u, </>) G H :— 
H\tt) x fl-V2(r) such that 

(aVu , V^)q + (2Uu + (£' - ±)0 , u) r = (/ , v) n + (fa + 2Uu , v) r , (3.3a) 

(ip , QJ0 — (& — i)«) r = -(^ , (£ - ±K)r, (3.3b) 

for all (v,i/j) G ft. 

The link between (3.1) and (3.3) is provided by u = u mt and </> = Vw cxt • v. Moreover, u ext 
is then given by the third Green's formula 

u ext (x) = &(u-u )(x) -V3<t>(x) forxGfi ext , (3.4) 

where the potentials 03 and R formally denote the operators 03 and .ft, but are now evaluated 
in f2 ext instead of V. Note carefully that we do not use a notational difference for the function 
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u G H 1 ^) and its trace u G H 1//2 (r), for which we compute the boundary integrals Wu and 
(£- \)u in (3.3). 



3.2. Existence and uniqueness of solutions. Assumption (3.2) guarantees pointwise 
ellipticity 

(v — w) for all v, w G 1R and x £ Q (3-5) 
with (•) denoting the Euclidean scalar product on IR d , as well as pointwise Lipschitz continuity 
\%(x)v - 2l(x)w| < C lip \v - w\ for all v, w G M d and x G fi, (3.6) 

1/2 

where C mon = C m i n , C up = Cmax > do not depend on x. 

Existence and uniqueness of the solution u = (u, 4>) of (3.3) rely on the J f/" 1 / 2 (r)-ellipticity 
of the simple-layer potential 33. Details are found e.g. in [CS2, AFFKMP]. In 3D, the simple- 
layer potential 2J is always elliptic, i.e. 

IMIir-v-cr) S ^ > *Wr far all ^ G F" 1 / 2 ^). (3.7) 

To ensure ellipticity in 2D, it suffices to scale flcK 2 appropriately so that diam(f2) < 1. 

In 2D, the radiation condition u(x) = 0(l/\x\) either has to be relaxed to u(x) = C(log \x\) 
as |x| — > oo or the given data must satisfy the compatibility condition J n fdx + J r 0o dT = 0. 
The latter implies J r 4>dT = and hence the right decay (3.1e) of w ext at infinity, as can be 
seen from (3.4), cf. e.g. [S, Section 6.6.1]. 

3.3. Galerkin discretization. Let Sf be a K-shape regular triangulation of the coupling 

boundary T in the sense of Section 2.1. In addition, let Ti be a regular triangulation of Q 

into compact and non-degenerate simplices T G %. As above, we assume that Q as well as 

?r 

't ' 

of 



T are exactly resolved by % and £ J, and K-shape regularity of % is understood in the sense 



diam(T) a! . . 

o\Tt) := max — < k < oo (3.8) 

with | • | denoting the usual volume measure on M. d . By Sp, we denote the set of facets of % 
which lie inside of Q, but not on the coupling boundary. 

For the discretization, we use conforming elements and approximate a by a continuous 
7rpiecewise polynomial Ug G S p {%) C H l {Vt) of degree p > 1. Moreover, (ft is approximated 
by a (possibly discontinuous) £j-piecewise polynomial G V q (£f) C iJ _1 / 2 (r) of degree 
g > 0. We stress that the usual link between p and q is g = p — 1, and the lowest-order case 
would be p = 1 and q = 0. 

The discrete spaces read 

X t := 5 P (7^) x C H\n) x iJ" 1 / 2 ^) = U, (3.9) 

where the product space H, equipped with the canonical norm 

IIMH = {\Hh m + ll^ll^ 1/2( r)) 1/2 for v := (v^) G H, (3.10) 
becomes a Hilbert space. 
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The Galerkin formulation of (3.3) reads as follows: Find U£ = (U£, G Xt such that 

(SlVEtf , W*)n + (WU* + (# - , r £ )r = (/ , V«) n + (0o + 2H«o , V«)r, (3.11a) 

, 9J$J - (£ - |)C/;) r = §)u )r, (3.11b) 

for all V^ = (V«,^)eAi. 

Again, it is known that (3.11) admits a unique discrete solution U£ G which is quasi- 
optimal in the sense of the Cea lemma 

|||u - U^lll < C C ea min III vi - V,|||, (3.12) 

where the constant Ccea > depends only on Q, see [AFFKMP] resp. [AFP, Appendix] for 
the fact that, contrary to [CS2, Corollary 3], no additional assumption on % or Sf is needed. 

3.4. Perturbed Galerkin discretization. The right-hand side of the discrete formula- 
tion (3.11) involves the evaluation of 2Uu and Ru , which can hardly be performed ana- 
lytically. Moreover, so-called fast methods for boundary integral operators usually deal with 
discrete functions, cf. [RS]. Therefore, we propose to approximate at least the given bound- 
ary data u G iJ 1 / 2 (r) by appropriate discrete functions and proceed analogously to [KOP]: 
To that end and to provide below a local measure for the approximation error, we assume 
additional regularity u G H 1 ^) and use the Scott-Zhang projection Ji : L 2 (T) — > S p (Sf) 
from [SZ] to discretize Uqj = JiUq G S p (Sf). Now, the perturbed Galerkin formulation 
reads as follows: Find = (Ui, <Ev) G Xi such that 

(WU e , VV e }n + (WBUt + (# - l)®t , V e )r = (f , V t )a + (0o + WU 0>t , V t )r, (3.13a) 

93$* - (£- §)^)r = -(^, (it- l)U 0/ ) T , (3.13b) 

for all \~£ = (Ve, ^) G X#. Compared to (3.11), the only difference is that (3.13) involves the 
approximate data Uq/ instead of no on the right-hand side. Consequently, the same argu- 
ments as before prove that (3.13) has a unique solution. Moreover, the benefit is that (3.13) 
only involves discrete boundary integral operators, i.e. matrices. 

Remark. We stress that the additional regularity assumption uq G H 1 ^) is also necessary 
for the well-posedness of the residual error estimator of Section 3.5. □ 



3.5. A posteriori error estimate. The overall goal of this section is to provide a residual 
a posteriori error estimate (3.18) for our discretization (3.13). To that end, we assume 
additional regularity no G H 1 ^) and <po G L 2 (T). Recall the notation %, Sf 1 , and Sf from 
Section 3.3. We define the volume residual 

Vi(Te) 2 = ^( T ) 2 ' where ^( T ) 2 = \ T \ 2M 11/ + div(2lW,)||i 2(T) , (3-14) 

TeT f 

the jump contributions across interior edges 

Vi(£?) 2 = J2vi(E) 2 , where m {Ef = |£|V(<*-i) || [(2lVf/,) • v\ ||| 2(£) , (3.15) 
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and the boundary contributions on the coupling boundary 

n^IT = E ^ E ) 2 ' where ^ 2 = ^(e? + v?\Ey, 



4 1] (e) 2 = iei 1 ^ ||0o - (avc/,) ■ 1/ + 2U(tv - ^) + (I - 11 



vf\E) 2 = |£r /(d - 1} ||V r (9J*< - (| - - 17<" 

Moreover, we define the data approximation term 



(3.16) 



osc 2 



Eosc,(E) 2 , where osc,(£) 2 = |^|i/C^-D || Vr ( nQ _ U , e )\\h(Ey ( 3 - 17 ) 



Altogether, we thus obtain the following a posteriori error bound, whose local contributions 
are used in Section 3.7 below to steer an adaptive mesh refinement. 

Proposition 13. The error is reliably estimated by 



in tt ni2 / n 2 / n2 ... / n2 ) Ve(rf forreTeUS^ 

re ^ u£ r ^(r) 2 + osc,(r) 2 /orrG^. 

(3.18) 

T/ie constant C rc \ > depends only on n-shape regularity ofTe and Sf . 

Sketch of proof . Instead of solving the non-perturbed Galerkin formulation (3.11) of the 
weak formulation (3.3), we solve the perturbed Galerkin formulation (3.13) in practice. 
Put differently, = (Ui, $^) G Xt is the Galerkin approximation of the unique solution 
U£ = («£, <pe) G % of the continuous perturbed formulation 

(2lVu £ , Vv) n + (Wu e + (# - , u) r = (/ , v) n + (<Po + %W 0>t , v) T , (3.19a) 
(i/j , - (it - i)u,) r = , (£ - |)^)r, (3.19b) 
for all v = (v, if)) G By stability, the error between u and is controlled by 

1 /2 

|||u - u^lll < ||«o - U ,e\\ H i/2 {r) < \\h/ V r (% - U 0/ )\\ L 2 {r) ~ osc e , (3.20) 

where the required approximation estimate for the Scott-Zhang projection can be found 
in [KOP, Theorem 3]. Here, ht G V°(Sf) denotes the local mesh size function of Section 2.1. 
Therefore, it only remains to prove 



uj < m = ( E ^ T ) 2 + E %(^) 



, V2 

2 4 ' 



This follows along the lines of the 2D proof presented in [CS2, Theorem 4]. The only 
remarkable difference is that in the 3D case the estimate 

||<2J$ £ - (| - &){U QA - Ui)\\ H y2 {T) < \\hy 2 V r (^ e - (| - St){U 04 - U t )) |U 2(r) 

does not follow from continuity arguments as in [CS2], but from a Poincare-type estimate 
provided by [CMS, Corollary 4.2]. □ 
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Remark, (i) For the scaling of the different contributions of Qe defined in (3.14)-(3.17) 7 
recall that diam(T) ~ | CZ -1 1 1 / cf for a d- dimensional volume element T E 7} and diam(-E') ~ 
\E\ 1 ^ d " 1 ' > for a (d — 1) -dimensional boundary element E £ Sf or interior facets E £ Sf 1 . 
Altogether, the volume contribution (3.14) is weighted by diam(T), while all other contribu- 
tions (3.15)-(3.17) are weighted by diam(£ , ) 1/ ' 2 . 

(ii) In [KOP, Theorem 3], the approximation estimate in (3.20) is proved for any H x l 2 -stable 
projection Jn onto S p (Sf). For p = 1, it is proved in [KPP, Theorem 6] that newest vertex 
bisection guarantees H 1 -stability and hence also H 1 ^ 2 -stability of the L? -orthogonal projec- 
tion ILj : L 2 (T) — > S 1 (Sf) on the 2D manifold Y. In the case p = 1, we may thus also use 
Uq.i = TlgUo to discretize the given Dirichlet data, and Proposition 13 still holds accordingly. 

(iii) In 2D, one may also use nodal interpolation Uq/ = hu$ to discretize the data, and (3.20) 
holds accordingly. Details are found in [AFP, Proposition 1]. □ 

3.6. Local mesh refinement. Let Ti be a sequence of triangulations of Q which is obtained 
from a given initial triangulation 7o by successive local mesh refinement. We require the 
following assumptions on the local mesh refinement of the volume mesh with ^-independent 
constants < n < oo and < q < 1: 

(Tl) The triangulations % of Q are regular in the sense of Ciarlet and K-shape regular. 

(T2) Each refined element T £ % is the disjoint union of its sons T' £ li+i- 

(T3) The sons T' £ %+i of a refined element T £ % satisfy \T'\ < q\T\. 

(T4) The sons E' £ £f +1 of some refined facet E £ £f satisfy \E'\<q\E\. 

In addition, let Sf be a sequence of triangulations of the coupling boundary T which is 
obtained from an initial triangulation Sq by successive local mesh refinement. For the local 
mesh refinement of the boundary mesh, we assume the following: 

(El) The triangulations £f of T are regular in the sense of Ciarlet and K-shape regular. 
(E2) Each refined element E £ Sf is the disjoint union of its sons E' £ £e+i. 
(E3) The sons E 1 £ £f +1 of a refined element E £ £f satisfy \E'\ < q\E\. 

From (T2) and (E2), it follows that the discrete spaces are nested Xi C Xi + \. 

All of these assumptions are satisfied for the usual mesh refinement strategies, see e.g. [V]. 
For instance, one may use newest vertex bisection for both Tt and £f. We refer to e.g. [St] 
for details on the latter algorithm, but remark that (T3)-(T4) and (E3) are then satisfied 
with q = 1/2. Moreover, we stress that 2D and 3D newest vertex bisection only leads to 
finitely many equivalence classes of elements T £ IJ^o 

In the experiments below, we let £f = %\r be the induced triangulation of the coupling 
boundary T and use newest vertex bisection to refine %. In 3D, Sf can then equivalently 
be obtained by use of 2D newest vertex bisection from Sf. We stress, however, that this 
coupling is not needed in theory. 

3.7. Adaptive algorithm and convergence. In this section, we consider the following 
common adaptive algorithm, which steers the local mesh refinement by use of the local 
contributions of Q 2 = rj 2 + osc 2 of Proposition 13. 

Algorithm 14. INPUT: Initial meshes (%,Sq) for i := 0, adaptivity parameter 9 £ (0, 1). 

(i) Compute discrete solution ILj £ X^. 

(ii) Compute refinement indicators q^{t) from (3.18) for all r £ % U Sf U Sf. 
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(iii) Determine a set Aig C % U Sf U Sf which satisfies the Dorfler marking criterion 

T£TenM e Ee£^nM e E&£fnM e 

(iv) Mark elements T e 7^ D M.t and facets E 6 (£|P U £j) D .M^ for refinement. 

(v) Generate (7^+1, &?/ refinement of at least all marked elements and facets. 

(vi) Increase counter £ 1— >■ £ + 1 ; and goto (i). 

OUTPUT: Sequence of error estimators (geje^N ond discrete solutions (Ueje^. D 

Since adaptive mesh refinement does not guarantee that the volume mesh Ti and the 
boundary mesh Sf become infinitely fine, convergence XJe — > u is a priori unclear as I — > 0. 
However, we employ the concept of estimator reduction from [AFLP] to prove qi — > as 
£ — y 0. Finally, convergence — > u as I — > 00 then follows from the reliability (3.18) of Qg. 
These observations are stated in the following theorem. 

Theorem 15. Algorithm 14 yields a sequence of meshes %, S(_, and approximations XJi with 
the following properties: 

(i) there are constants < p < 1 and C2 > such that the overall error estimator q\ = 
rjf + osc 2 satisfies 

Q 2 t+l <p$ + C 2 (|||U m - U,||| 2 + ||^+i - U^f Hl/2{r) ) for all I > 0. (3.22) 

The constant < p < 1 depends only on the mesh size reduction constant < q < 1 
from assumptions (T3)-(T4) and (E3) and the adaptivity parameter < 6 < 1. The 
constant Ci > depends additionally on Q, n-shape regularity of Ti and Sf , and the 
polynomial degrees p, q. 

(ii) The error indicators as well as the oscillation terms tend to zero: 

lim rje = = lim osq. (3.23) 

l— >oo I— >oo 

(iii) The adaptive coupling algorithm converges: lim |||u — U^||| = 0. 

t— >oo 

Proof of Estimate (3.22) of Theorem 15. We recall that (a + bf < (1 + 5)a 2 + (1 + cT 1 )^ 
for all a, 6 £ R and arbitrary 5 > 0. 

• First, for the volume contributions (3.14) and the interior jumps (3.15), we argue as 
in [CKNS] . By use of the triangle inequality, we see 

r) £+1 (% +1 ) 2 = Yl |T'| 2/d ||/ + div(2lVf/ m )||i 2(T , ) 

T'eT £+1 

<(l + 5) \TT' d \\f + div(WU e )\\l 2{Tf) {32A) 

T'eT e+1 

+ (l + r 1 ) iTl^Hdiv^V^x-V^Hi^,) 

r'er f+ i 

The product rule gives div2t(VC/^+i - V^) = (divSt) • VU i+1 + 21 : D 2 (VE/ m - VU t ) with 
-D 2 (-) being the Hessian matrix. Using the triangle inequality and a scaling argument, the 
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arguments of the second sum in (3.24) are estimated by 
|T'| 2 / d ||div2l(Vf/ m -V^)||i 2(T0 

< |Tf / d (||div2t||i 00(T , ) ||V(t/ m - U,)\\\ 2[TI) + ||2l|||o O(T0 || J D 2 (VC/ m - V^)||i 2(T0 

< 1 1 21 1 1 ^.^(T') l|V(L^+i - Ut)\[\ 2 ^ T ,y 

The first sum in (3.24) is split into non-refined and refined elements, and assumption (T3) is 
used for the refined elements. Recall that assumption (T2) states that each element T G 7f 
is the union of its sons T' e Te+i- This gives 

E \T'\ 2/d \\f + d^vu e )\\l 2{Tl) 

< E \T\ 2/d \\f + dW(QlWU e )\\l 2{T) + q ^ d E |T| 2/d ||/ + div(2tVf/,)||i 2(r) 

TeT e nT e+1 TeT e \T e+1 

= m(%) 2 - (1 - q 2/d ) Tu(Ti\Ti+i)*. 
Altogether, this gives 

m + i{% +l f < (1 + 5)( Ve (T e f - (1 - q 2 ' d ) m{%\% +1 ) 2 ) 
+ {l + 8- l )C\\V{U l+1 -U i )\\l. m i 

where C > depends on Q, on /t-shape regularity of and the (local) W 1,00 -norm of 21. 
• Second, for the interior jumps (3.15), we argue 

m + i{sf +1 f= E \Ef^[mu e+1 -u]\\i 2{E/) 

<(i + s) E \Ef^\\[mu e -u)\\l 2m (3 _ 26) 



i + r 1 ) E l^l^llfstv^-^-HII^ 



E'P£p 



The arguments of the second sum in (3.26) are estimated by use of a scaling argument, 
namely, 

\Ef^\\[QlV(Ue +1 - Ut) ■ v]\\ 2 L . iEI) < ||2l||^.oc ( ^ +1B/) ||V(C/ m - U^^^ 

with w^+i^' = T| U f being the patch of E' — T' + r\T'_ e To treat the first sum 

in (3.26), the essential observation is that 21VL^ is continuous across new facets £g^-i\£g^ so 
that the respective jumps vanish. This and assumption (T4) for refined interior facets give 

E i^i^ii^-^iii,^ 

< |S| 1 /( d - 1 )||[2tV^-Hlli 2(B ) + 9 1/(d " 1) E l^| 1/(d_1) ll[2lVC/,-H 

= ^(^) 2 -(l-g 1 /(^i) ) ^(^^ +i)2 _ 
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Altogether, we obtain 

%i(d) 2 < (i + s)( m (s?f - (i - g^Ms^f) 
+ (i + r 1 )C||v(^ +1 -^)||i 2(n) , 

where C > depends on K-shape regularity of Ti+i and the (local) W /1,00 -norm of 21. 

• Third, we consider the first boundary contribution (3.16) of the discretization error 
estimator. By use of the triangle inequality, we see 

||0o - (2lW m ) • v + 2]J(tfy + i - U t+1 ) + (i - #)^ + i|| L 2 (E0 
< ||0o - lpXU t ) ■ v + 2U(fV - ^) + (| - #)^IU^') 
+ ||2lV(f/ m - U t ) ■ u\\v (JSI) + ||2U((f/ ,, +1 - U 0/ ) - (U e+1 - U e ))\\ L2{EI) 

+ r ll^+i - ^IU 2 (s') + ||#'(^+i - 
We sum these terms over all elements E' G Sf +1 and use assumption (E3) to see 
lEf/^ ||0o - (fXVUt) ■ v + W(U 0/ - U t ) + (| - i^llW) 



< ^ |^| 1 /(^i)||0 o _ (avc^) • v + 2S(cv - u t ) + (f - 



g i/(d-D ^ i^i 1 /^!) || 0O _ (av^) ■ + mj(u , e - m) + (| - &)$ e 



'£\\l 2 (E)- 

Ee£f\£f +1 



A scaling argument reveals 

Dllofw/rr tt \ _ . 1 1 2 llvr/Tr rr M, 

lL 2 (f7); 



£ |E'| 1/(d - 1) ||2tV(t/ m - [/,) • v\\ 2 L2{EI) < \\V(U £+1 - U e )\tf 



where the implied constant depends on the (local) W 1,00 -norm of 21 and K-shape regularity 
of 7^+i . The inverse estimates of Corollary 2 for 8! and 233 shows 

\Ef/^(\\w((u 0/+1 - u 0ti ) - (u e+1 - u e ))\\i HEI) + \\#(<s> £+1 - mi^)) 



E'e£f +1 



= \\h\^<m({u w - u 0/ ) - (u e+1 - u t )) + H^O'(^ + i - ^)||| 2(r) 

< ||(^+1 - U 0>i ) - {U £+1 - ^)||^ 1/2(r) + - $£||^-i /2(r) . 

An inverse estimate from [GHS, Theorem 3.6] proves 

£ \Ef/( d -V\\$t +1 - ^\\ L , m = - m^r) S H $ m - *i\\ 2 h- 



E'e£f +1 



i/2(r)) 
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where the implied constant depends only on T, the K-shape regularity of £f +1 , and the 
polynomial degree q. Combining all these estimates, we obtain 

4U^U 2 < (1 + *)(tfW) 2 " (1 " ? 1/(d - 1) )% (1) (^ r + i) 2 ) 

+ (1 + 5^)C (\\V(U t+ i - U e )\\l 2{ ^ + \\U e+1 - U t f Hl/2{v) (3.28) 

+ H^W+l — U 0> £ ||#l/2( r ) + \\$t+l — ^ll/f-l/2(r)) 

for the first contribution to f]e{£f) from (3.16). The constant C > depends only on T, on 
K-shape regularity of £f +1 and Te+i, and on the polynomial degree q. 

• Fourth, we consider the second boundary contribution (3.16) of the discretization error 
estimator. Arguing as before, we now use the inverse estimates from Corollary 2 for QJ and 
R as well as a local inverse estimate 1 1 /i^Vr 1 1 z, 2 (r) ^ HV^Iifi/sfr) f° r an S p (%) which is 
e.g. found in [CP2, Proposition 4.1] for 2D and [AKP, Proposition 3] for 3D. We obtain 

rifUeU 2 < (1 + 8)(v?>(£fi 2 - (1 " <l lKd - 1) )vt\^Uf) 

+ (1 + 5~ 1 )C (\\U t+1 - U e \\ 2 Hl/2{r) (3.29) 

+ H^W+i — U 0) i\\ 2 Hl/ 2^ + \\$t+i — ^ll/f-i/2(r)) 

for the first contribution to rje(£f) from (3.16). The constant C > depends only on K-shape 
regularity of £f +1 and Te+i, on T, and on the polynomial degree p. 

• Fifth, we consider the oscillation terms from (3.17). Arguing as in the previous steps, it 
follows that 

osc 2 +1 <(l + 5) (osc 2 - (1 - q^)osc e (£f\£f +1 ) 2 ) 

+ (1 + 8 1 )C \\Uo t £ + i — Uo,i\\jii/2(r\, 

where the constant C > depends on T, K-shape regularity of £f, and the polynomial degree 
p. 

• Sixth, we combine the reduction estimates (3.25), (3.27), (3.28), (3.29), (3.30) for the 
different parts of Qg. By a norm equivalence (with constants depending only on Q) we have 

l 2 (n) + ll^m — Ue\\ 2 Hl/2 ^ — \\Ui+i — ^11^1(0)- Together with < q < 1 and 
hence q 2 ^ d < q 1 /( d ~ 1 \ we obtain 

Qe+i = m+i(%+i) 2 + r) e +i(£?+i) 2 + Ve+i(£e+i) 2 + °sc 2 +i 

<(l + 5) (^(T,) 2 + m {£?) 2 + Ve (££) 2 + osc 2 

- (i - q 1/id - 1] ) H%\T l+1 f + m{£f\£f +1 ) 2 + m(£l\£f +1 ) 2 + osc^AO^ 



+ (1 + 5 1 )C (||£^+i - Ue\\ 2 H i^ + \\U 0) i + i - U 0) £\\ 2 H i/2^ + \\$t+i - ^lln-va^))' 

and the constant C > depends only on Q, K-shape regularity of Tg and £f, and the 
polynomial degrees p, q. To abbreviate notation, we define the index set Zg := % U £f U £f. 
Together with norm equivalence \\U e+1 - U e \\ 2 Hl ^ + \\$t+i - ^e\\ 2 H -i/2r T ) — |||U^+i - U e \\\ 2 , 
the previous estimate takes the form 

ej +1 < (1 + 5){g 2 - (1 - q 1 ^) 6i (Z e \Z i+1 ) 2 ) 

+ (1 + 5^)C (|||U m - U,||| 2 + \\U 0/+1 - U 0/ \\ 2 Hl/2{r) ). 
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• Seventh, recall that marked elements and facets are refined so that Aii C Xi\Li + i. With 
the abbreviate notation, the Dorfler marking (3.21) takes the form 

eg] < Qe(Me) 2 < g e (l e \l e+1 ) 2 . 

Therefore, the estimate (3.31) becomes 

Q 2 l+l = (1 + 6)(1 - (1 - gVC^i))^ Q \ + (i + s-^C (|||U m - U,||| 2 + \\U 04+l - U 0/ \\ 2 Hl/2{r) ). 

Now, note < (1 - q l/{d ~ l) )6 < 1. Therefore, we may choose some sufficiently small 
5 > such that p := (1 + 5)(l - (1 - g 1 /^" 1 ')^) < 1. This concludes the proof of the 
estimate (3.22). □ 

We recall that assumptions (T2) and (E2) imply X^ C Xg + i. The proof of conver- 
gence (3.23) of the adaptive coupling relies on the following a priori convergence result 
from [AFP, Proof of Proposition 10] which essentially follows from the nestedness Xi C Xi + i 
of discrete spaces and the validity of Cea's lemma (3.12). 

Lemma 16. According to assumption (T2) and (E2), the limit := lim^ U£ of Galerkin 
solutions exists strongly in 7-L. If the limit Uq >00 := lim^ Uoj exists strongly in if 1//2 (r), then 
also the limit Uqo := lim^U^ of the perturbed Galerkin solutions exists strongly in 7-L. □ 

Moreover, the following result from [FPP, Lemma 18] proves a priori convergence of Uq^ 
for the Scott-Zhang projection. We remark that [FPP, lemma 18] is stated and proved for 
p — 1 and if 1 (f2), but the proof holds without any modification also for general p > 1 and 
H 1 ^). We recall that assumption (E2) implies nestedness S p (£f) C S p (£f +1 ), whereas (El) 
states uniform K-shape regularity which enters the approximation and stability estimates of 
the Scott-Zhang projection. 

Lemma 17. According to (El) and (E2), the limit J^g = \im.i Jig of the Scott-Zhang 
projections exists strongly in H 1 ^) for all g G H l {T). □ 

Proof of Convergence (3.23) of Theorem 15. According to Lemma 17, the limit £/o,oo = lim^ Uo t e 
exists strongly in H l (T) and hence also in i7 1//2 (r). Therefore, Lemma 16 applies and proves 
that Uqo := lim^U^ exists in %. Consequently, the estimator reduction estimate (3.22) may 
be written as 

Qt+i < P Qi + oli for all € > 

with a non-negative zero sequence on — » as i — > oo. It is thus a consequence of elementary 
calculus that lim^ qi = 0, cf. [AFLP, Lemma 2.3]. With q 2 = rjj + osc 2 and reliability 
lll u — U^IH < Qe, we conclude the proof. □ 



Remark, (i) For p = 1 and 2D, one may also use nodal interpolation Uq i = I#uq to dis- 
cretize the Dirichlet data. It is proved in [AFGKMP, Proof of Proposition 5.2] that the limit 
Uq.oo = lim^f/o^ then exists strongly in H 1 ^) and hence also in H l l 2 (Y). 

(ii) For p = 1 in 3D and newest vertex bisection, one may also use the L 2 ' -orthogonal pro- 
jection Hi : L 2 (T) — > S 1 (Sf). It is proved in [KPP, Theorem 3] that IL is H 1 -stable and 
in [KOP, Lemma 12] that the limit Uq j00 = lim^ Uqj therefore exists weakly in H 1 ^) and, by 
interpolation, hence strongly in if 1//2 (r). 

(iii) In either case, the claims (3.22) -(3.23) of Theorem 15 remain valid. □ 
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3.8. Extension to nonlinear transmission problems. We consider the interface prob- 
lem (3.1) with nonlinear 21 : R d -> R d and 2lVtt(x) := 2l(Vit(x)), see e.g. [CS2, AFFKMP]. 
Under monotonicity (3.5) and Lipschitz continuity (3.6), the coupling formulation (3.3) as 
well as the discretizations (3.11) and (3.13) still admit unique solutions. Moreover, the 
a posteriori error analysis of Proposition 13 remains valid. Here, note that 21VL^ is Tit 
piecewise Lipschitz continuous so that all occurring terms as e.g. the jump terms (3.15) are 
well-defined. 

While the convergence result of Theorem 15 holds for arbitrary polynomial degrees p > 1 
and q > in case of linear problems, we did not succeed to prove the same for nonlinear 
transmission problem. Difficulties are only met in the first step of the proof. In (3.24), the 
term 

|T / | 2/d ||div(2lV^ +1 - 2tW,)||i 2(T , } ~ \\h e div e (mU e+1 - 2lV^)||£ 2(n) 

arises, with div^ the Trpiecewise divergence operator. However, for lowest-order volume 
elements p = 1, this term vanishes. Thus, our proof also covers the adaptive FEM-BEM 
coupling for certain nonlinear coupling problems and lowest-order elements p = 1, where 
21 : R d ->• R d satisfies (3.5)-(3.6). 

For p > 1 and 21 being linear, we used a scaling argument to bound the latter term in (3.24) 
by ||2l|| || Vt/^+i — VL^||i2(n). For 21 : R d — > R d being nonlinear, we did neither succeed 
to find a similar bound (which seems, in general, to be impossible) nor to prove that this 
term vanishes as £ — > oo. 

4. Numerical experiments 
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Figure 1. Z-shaped domain and initial triangulations 7o, £q with #7o = 14 
triangles and =10 boundary elements. 

In this section, we present 2D calculations for the perturbed Galerkin discretization (3.13) 
of the model problem (3.1). We consider different choices of the polynomial degrees p and q 
for the spaces S p (7~i), V q (£f). For ease of implementation and also for stability reasons, we 
restrict ourselves to the case E\ = 7e\r, i- e -> the boundary mesh is taken as the trace of the 

27 




10° 




10 number °of volume elements 1 V 1 ° number of boundary elements M 

Figure 2. err £ (fi) and Qz{Vl) vs. N = (left) as well as err^(r) and g e (T) 
vs. M = #£f for adaptive and uniform mesh refinement with p = 2 and 5=1. 
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Figure 3. err^(fi) and Qi(ti) vs. N = #Te (left) as well as err^(r) and ^(r) 
vs. M = #£f for adaptive and uniform mesh refinement with p = 2 and q — 0. 



volume triangulation. All computations were performed on a 64-BIT Intel(R) Core(TM) i7- 
3930K Linux work station with 32GB of RAM. For the computation of the discrete boundary 
integral operators we used the MATLAB BEM-library HILBERT [ABEM], and all systems of 
linear equations were solved with the Matlab backslash operator. 

We employ the adaptive Algorithm 14 with 9 = 0.25 and compare the results with uniform 
mesh refinement (this can be realized by setting 9 = 1 in Algorithm 14). The domain O is 
taken as the Z-shaped domain visualized in Fig. 1. We take 21 = Id in (3.1) and prescribe 
the data /, u , O such that the exact solution is given by 

v} nt (r,ip) =r 4/7 sin( fv?), 

ext. \ x + y + 0.25 
u (x, y) = ,u 2 ,, , !, 2 ; 
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10 number of volume "elements 1 V 1 ° number of boundary elements M 

Figure 4. err £ (fi) and ge(Q) vs. iV = j^Ti (left) as well as err^(r) and ^(r) 
vs. M = #£f for adaptive and uniform mesh refinement with p — 1 and 9 = 1. 




Figure 5. err e (Q) and gt(tt) vs. A 7 = #Te (left) as well as err^(r) and g e (T) 
vs. M = #£f for adaptive and uniform mesh refinement with p = 1 and q = 0. 



the polar coordinates (r, (p) are taken with respect to the origin (0,0). The prescribed 
solution u mt has the typical singularity at the reentrant corner (x, y) = (0, 0), which leads to 
a reduced order of convergence (9(/i 4//r ) for uniform mesh refinement. 

Recall that the a posteriori error estimator qi from (3.18) is split into volume contribu- 
tions (3.14)-(3.15) and boundary contributions (3.16)-(3.17) 

Q \ = Qi (n) 2 + Ql {Yf. 

Arguing as in [AFP], the Cea-type quasi-optimality produces 

|||u - Ui\f = \\u - U e \\ 2 m{n) + ||0 - ^||^-i/ 2(r) < ||« - ^||hi(q) + Wh 1 / 2 ^ - ^)||£ 2(r) 

=: (err,(fi)) 2 +(err,(r)) 2 

and thus provides a computable upper bound for the overall error. We plot the error contri- 
bution en>(f2) as well as the corresponding estimator part Qe(Q) versus the number of volume 
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elements N = jfH~i- We also plot the boundary contributions err^(r) and ge(T) versus the 
number of boundary elements M = #£f. We observe convergence rates proportional to 
N~ a resp. M _/3 with some a, (3 > for all computed quantities. The optimal convergence 
rate for the FEM part with S p (7~e) is a = p/2, whereas the optimal rate for the BEM part 
with V q (£f) is (3 = 3/2 + q. Consequently, the optimal overall convergence rate for the 
FEM-BEM coupling is proportional to N~ a + M _/3 . 

The use of a uniform mesh refinement leads to suboptimal convergence rates for both the 
volume and boundary quantities. In particular, we observe the rate a = 2/7 independently of 
the chosen polynomial order p, q, see Figures 2, 3, 4, and 5. We recall that for uniform mesh 
refinement there holds N ~ M 2 and therefore the overall convergence rate is proportional 
to N- 2 ' 7 . 

In contrast to the uniform approach, the adaptive mesh refinement strategy recovers the 
optimal overall convergence rate. We observe optimal rates a = p/2 and (3 = 3/2 + q for 
p = 2, q = 1 as well as p = 2, q = and p = 1, q — 0, see Figures 2, 3, and 5. In the 
case of p = l,q — 1 we observe that the error estimator contribution Qe(T) converges with 
order M _3//2 , whereas the error quantity err^(r) has some slightly higher convergence rate. 
Nevertheless, we stress that the optimal overall convergence order is achieved, since our 
computations also show that N ~ M 2 in the case of p — 1, q — 1 and therefore the overall 
error is dominated by the FEM part, which converges with order a = 1/2, see Figure 4. 

In conclusion, our numerical experiments underline the fact that the adaptive mesh re- 
finement strategy of Algorithm 14 achieves the optimal convergence rates in the presence of 
singularities of the exact solution. 
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